Contribution of the screened self-energy to the Lamb shift of quasidegenerate states 
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Expressions for the effective Quantum Electrodynamics (QED) Hamiltonian due to self-energy 
screening (self-energy correction to the electron-electron interaction) are presented. We use the 
method of the two-time Green's function, which handles quasidegenerate atomic states. From these 
expression one can evaluate energy corrections to, e.g., ls2p 3 Pi and ls2p 1 P\ in helium and two- 
Q , electron ions, to all orders in Za. 



In the last ten years, experiments in the spectroscopy of helium [Jl] have become two orders of magnitude more 
precise than the best theoretical energy level calculations available (see, e.g., Refs. [|]j^]and references therein). Several 
experiments are now focusing on Helium and hcliumlike ions ls2p 3 Pj fine structure |§-|l§], with the aim of providing 
a new determination of the fine structure constant and of checking higher-order effects in the calculations. In this 
case the theory is again a limiting factor. In this context a direct determination of all a 2 contributions to all order 
in Za is necessary to improve reliability and accuracy of theoretical calculations (a being the fine structure constant, 
O ' an d Z the charge of the nucleus). 

A difficulty in the study of the (ls2p 1 / 2 )i and (ls2p 3 / 2 )i levels is that they are quasidegenerate for low and middle 
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Z ions |13]; this precludes the use of the Gell-Man-Low and Sucher method [14 15 1 to evaluate QED energy shifts of 



atomic levels. In fact, this method has two important drawbacks: it does not handle quasidegenerate energy levels, 
and it leads to a difficult renormalization procedure when applied to degenerate states. (The latter problem has only 
been tackled up to second-order in a Jl6| , p^[ .) 

We use the method of the two-time Green's function 18 -20), rigorously derived from QED (for the most detailed 
description of this method, see ) . To the best of our knowledge, only the method recently proposed by Lindgren [^2| , 
closely modeled to multireference-state Many-body perturbation techniques, is designed to work for quasidegenerate 
' states. 

We evaluate the contribution of the screened self-energy diagrams 
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to quasidegenerate energy levels in heliumlike ions. Our results can be easily extended to ions with more than two 
\ electrons along lines similar to those found in p^] . 

First approximate evaluations of the contribution of these diagrams for isolated states in two- and three-electron 
ions were performed in Refs. 24- 27|. Accurate calculations from the first principles of QED were accomplished in 



Refs. ]2q-|3(| for the ground state of heliumlike ions and in Refs. |31 32 for the 2s and 2pi/2 states of lithiumlike ions. 
The other two a 2 corrections to the electron-electron interaction have also been calculated for isolated states in two- 
and three-electron ions: the vacuum-polarization screening J13|J29|pO|,|3^;p4T| , and the two-photon exchange diagrams 
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[p5| |38[ . In [|13[ , the vacuum polarization screening for quasidegenerate states of heliumlike ions was evaluated as 
well. Some results for the direct contribution of the self-energy correction to the Coulomb interaction are also 
available @||. 

As depicted in diagrams ([!]) , the interaction between the two electrons through photons is treated perturbatively. On 
the contrary, the binding to the nucleus is included non-perturhatively in the method we use, since the corresponding 
coupling constant is Za. Such a treatment is obviously mandatory for highly-charged ions. Furthermore, it allows 
one to compare non-perturbative (in Za) results to (semi-)analytic expansions in Za (see [ pt)[ for a review). 

We derive the effective (finite-sized) matrix hamiltonian H , whose eigenvalues give the contribution of QED to a 
group of energy levels [ p3[ . The diagonal entries of the hamiltonian that we evaluate correctly reproduce previous 
expressions of the screened self-energy, while the new, non-diagonal entries that we derive allow one to obtain a 
second-order QED correction to quasidegenerate or degenerate energy levels. 
Relativistic units h = c = 1 are used throughout this paper. 

If we have s quasidegenerate energy levels the effective hamiltonian H is an s x s matrix restricted to these 

levels p3fl. Let us introduce some notations in order to express this hamiltonian. The second-order contribution 
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to this hamiltonian H = if W + ffW + + ... is constructed from a projection matrix P and an energy matrix 

H m = K W _ I { p(D 5 ^(i) } _ I {P (2) 7 ^(o) } 

+ |{[P«]2 j A-(0)} + IpWkWpW, (2) 
8 4 

where the notation {,} represents the usual anticommutator, and where the superscripts indicate the number of 
photons of the diagrams that contribute to each term of the perturbative expansion P = P'°) + PW + . . . and 
K = + + . . .; the s x s matrices P and K, which are defined as |2C|] : 

P=-LL Eg {E) (3a) 
2m J p 

K=^-JdEEg(E), (3b) 

where g{E) is the sxs matrix restriction of the Green's function to the s unperturbed atomic levels under consideration, 
and where T is a contour that encloses each of the Dirac atomic energy levels with a positive orientation p3j . 

We directly evaluate the hamiltonian matrix elements of Eq. (^) between states of different energies 
and put them in a form that readily displays the limiting case of identical energies; we checked by a direct calculation 

(2) 

of the diagonal matrix elements that they can be obtained from non-diagonal elements H ml , by taking the formal 
limit En^ — » . All the subsequent derivations of -ffl?/ will thus be done with E^ ^ E^/ . 

The first diagram of (|l|) appears only in the second-order matrices and P^ in Eq. (j|). As usual, we must 
calculate a reducible and an irreducible contribution; as can be seen in subsequent calculations, it turns out that 
the correct extension of these notions to quasidegenerate states is the following: in the first diagram of Eq. (Q), 
the contribution of intermediate electrons with a Dirac energy Ek such that Sk + s n ' coincides with one of the s 

P (2) 

energy levels under consideration and must be separated out from the contribution of the other intermediate electron 
states; the first contribution (called reducible) requires a different mathematical treatment from that of the second 
contribution (called irreducible). 

Thus, the irreducible contribution is obtained by summing over almost all electron states k in the first diagram of 
Eq. (0); we first show that it is sufficient to remove only one state k from the sum over states in the first diagram of Eq. 
(p]). We see that an intermediate energy Ek + £ n > t can coincide with an unperturbed atomic levels E^ s only if the 

R (2) 

electron k has the same principal quantum number as the electron n' p ,,y, on the other side of the self-energy, because 
otherwise the total energy Ek + £„' would lie largely out of the range spanned by the unperturbed quasidegenerate 

P (2) 

energy levels located around E^P) = e r j + £„< 

n P'(l) -P'(2) 

There is an additional selection on the electrons k to be removed: since the total angular momentum, its projection, 
and parity are conserved by the self-energy operator E [Eq. (^) below] , as can be seen by integrating over angles using 
standard techniques Ej ], the contribution of electrons k that do not share the same quantum numbers (k, m) as the 
electron n' p ,,^ in the first diagram of Eq. ([!]) is exactly zero. 

We denote the individual electrons of a state n by n\ and ri2, in an order which is arbitrary but that must remain 
fixed. With these notations, our evaluation of the irreducible part of the first diagram of ([!]) to the effective hamiltonian 
(0) takes a simple form and reads (Dirac energies are still denoted by £&): 



rjscr. SE, irr. 
H nn> 



(" 1 ) PP '( ( n P(i)\Z{e np(1) )\k) — — (kn P( 2)\I(£ npm -e n > )\m[m' 2 ) 

p,p> W^m ™ p(1) k 



+ £ (n P (i)n P(2) \I(e npw - £n' p , (1) )\n' P , (1) n' pl{2) )- — — (fc|£(e n , p/(i) )|n^ (1) ) J 

+0[a\E^-E^% (4) 

where (— l) pp is the signature of the permutation P o P' (P and P' arc permutations of {1,2}.), where the sum 
over k is over (almost) all possible intermediate Dirac states, and where the photon exchange and the self-energy of 
diagrams (P are represented by the following usual operators p2j : 
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(ab\I(u)\cd) = e 2 J d 3 Xl J d 3 x 2 [^(^1)0^(^1)] 

x[i(;l(x2)a u ip d (x2)}D^(uj;x 1 - x 2 ) (5) 
k\I(u)\kb) 

io)-(p-uy 



, ,™ , I j i 1 f , sr^ (ak\I(uj)\kb) 



in which e is the charge of the electron, a M = (1, a) are the Dirac matrices, and where ip denotes a Dirac spinor; the 
photon propagator D is given in the Feynman gauge by 



exp 



(i\r\^fuj 2 - p 2 + io) 



D vv ,{w;r) = g vv > ^ — — (7) 

47r|r| 

where fi is a small photon mass that eventually tends to zero, and where the square root branch is chosen such as to 
yield a decreasing exponential for large real- valued energies lu. 

The last term in Eq. (Q) represents a contribution of order a 2 which is multiplied by a factor that tends to zero as 

can be shown (see Ref. that such a term does not contribute to order a 2 and that it can 

therefore be omitted. 

We note that result (jZ]) readily yields diagonal elements by taking the (formal) limit E n ^ — E^) — > 0. 
The hamiltonian (|J) contains the contribution of many first-order diagrams through the operators P^ 1 " 1 and . 
We must consider here the contribution of the photon exchange and of the self-energy 
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n P'(2) 



their contribution to Eq. (g) cancels a part of the reducible screened self-energy. We thus evaluate in the following 
the contribution of both diagrams of Eq. (g) to the terms -\{P {1 \ K^} + §{[P (1) ] 2 , K (0 *>} + |p(l) (o)p(l) f the 
effective hamiltonian. 

The energy and projection matrices K and P of Eq. ([j|) have been calculated for the photon-exchange diagram 
in Jl3| ] ; this allows one to evaluate any integral due to the photon exchange that appears in the effective hamiltonian 

In order to derive the contribution of the one-electron self-energy, let us show that the evaluation of the self-energy 
contributions to the hamiltonian (||) boils down to the calculation of contour integrals of the form 

l dEg s n f(E) and -L J d E E g™(E) (9) 

where gf^(E) are diagonal elements of the self-energy Green's function; in other words, the contour T that surrounds 

all the levels in Eq. (g) can be replaced by the contour that surrounds En only, and non- diagonal elements of the 
self-energy Green's function are not relevant. The contour integrals of Eq. (]|) have both been evaluated in |p2| , so 
that no further quantity is required in order to obtain the self-energy contribution to the hamiltonian (^|). 

Let us prove the above statements. As mentioned before, angular momentum conservations constrain the self-energy 
operator E to be zero between states with different angular quantum numbers (k, to); and since the atomic levels 
we consider have the same principal quantum number (they are quasidegenerate), the self-energy Green's matrix is 
diagonal: 

gf*,(E) = iin^n', (10) 

where n and n' are the sets of quantum numbers of two of the s levels under consideration. 

Furthermore, the Green's function g n ®(E) has only one pole inside the integration contour T, namely at E = E n °\ 
Therefore, integrating over the full contour T in the hamiltonian (|^) amounts to integrate over the contour T n that 

surrounds only E n , since the Green's function is analytic inside the contours that encircle the other energies. 

We thus see that the contribution of the self-energy to Eq. (0) depends only on contour integrals of the form (||), 
which are known analytically p3| . 

With the help of some published analytical formulas, we obtain the following contribution of the photon exchange 
(see Eqs. (27) and (28) in and of the self-energy (see Eqs. (36) and (37) in [3^]) to the effective hamiltonian (||): 
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where S' represents the derivative of the self-energy operator (^J) with respect to the energy that flows in it, and 
where the two possible energies for the photon in the photon-exchange diagram are Ai = e 
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e n i and A2 

-p'(i) 



'™P(2) 



As seen above, the reducible part of the first diagram of Eq. (|l|) represents the contribution of an intermediate 
(For the second diagram, the reducible part is similarly obtained through an intermediate 



electron k 



l p'(i)- 



electron k = n P ^ .) The evaluation of the reducible contribution follows steps similar to those used for the irreducible 
part. The contribution of diagrams (^J) to the effective hamiltonian H&> , which is g iven in Eq. (pi]), cancels a few terms 
of the contribution of the reducible diagram, as for diagonal matrix elements |32j| ; the total reducible contribution to 
Eq. (0) is then found to be: 
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(12) 

where d x \ Xq represents the derivative with respect to x at the point xq. 

For the vertex diagram [second diagram of (EI)], the two-time Green's function method yields the following contri- 
bution to (||): 
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with the same notations as before; the sum is over all pairs of Dirac states. 

We thus have obtained the full contribution [Eq. (|J) + Eq. (|l2])+Eq. (|l3|)] of the screened self-energy diagrams (jl|) 
to a finite-sized effective hamiltonian which acts on a few atomic energy levels (in the general case: quasidegenerate, 
fully degenerate or isolated); the eigenvalues of this hamiltonian give the QED prediction for the energy levels. We 
have also taken into account the contribution of the first-order diagrams (||) to the second-order hamiltonian (|J). 

The results presented here extend previous derivations of the screened self-energy contribution to the Lamb shift, 
which were restricted to the evaluation of the energy shift of an isolated level. The diagonal terms of the effective 
hamiltonian that we have evaluated confirm previously published results. The new, non-diagonal matrix elements 
of the hamiltonian that we obtained allow one to calculate the energy shifts of quasidegenerate levels and to extend 
numerical calculations [^4|,^8 - 3lJ 4^] to such levels. 
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